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Introduction 


The Mariner 9 orbiter spacecraft observed that the 
Martian atmosphere was unexpectedly warm during the long- 
lived global dust storm of 1971 o Subsequently, Gierasch 
and Goody (1972) noted that this was consistent with the 
direct absorption by airborne dust of a significant fraction 
of the incoming solar radiation. This suggested a simplistic 
radiative heating model in which the heating per unit mass, 

0^, is given by pure absorption of some fraction, of the 
direct solar beam, Conrath (1975) has used such a model to 
examine the dissipation of the Martian 1971 -dust storm. 

If the dust particles are assximed to be spheres, Mie 
theory can be used to include the effects of radiative scat- 
tering by the dust if the size distribution arid the wavelength 
dependent refractive index — particularly its complex component 
are known. Moriyama (1974, 1975) has computed the solar heat- 
ing and infrared cooling of the Martian atmosphere for various 
degrees of dustiness by assuming that the Martian dust parti- 
cles are quartz particles having the same size distribution 
as terrestrial dust observed over northwest India. 

Recently, however, analysis of the Mariner 9 ultraviolet 
spectrometer (UVS) observations have yielded new constraints 



2 


on the size distribution and on the complex refractive index 
of the airborne Martian dust. In this paper we estimate from 
various sources the complex refractive index in the solar 
spectral region (0.2 < X < 5.0pm). With this index and the 
size distribution estimated from the UVS data (Pang ^ , 

1976) we use Mie theory to compute the radiative parameters 
requirejd- fox. ineliiding. radietive. scattering, in simplified 
radiative transfer approximations. These approximations 
then enable us to compute the wavelength-integrated solar 
heating for various dust- laden Martian atmospheres. 

We have deferred to an appendix material discussing the 
relative merits of the delta-Eddington (Joseph ^ al . , 1976), 
four-stream discrete ordinate (Liou, 1974) and four- component 
multiple stream methods for the optical parameter range ap- 
propriate to the Mars global dust storm. The calculations 
presented in the main body of the paper were computed using 
the delta-Eddington approximation . 
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II, Complex -Refxactive Index for Solar Wavelengths 

T-o include scattering in the radiative transfer equations, 
we need to know Che normalized phase function F(cos ©) which 
describes the probability that a photon will be diverted by an 
angle 0 from its original path during a single encounter with a 
nona-b sorbing d-u-&t- part-icJLe.. Eox s c a.t.t,erL±n^- which, is azimuthally 
independent, P(.cos S) can he written where ^ and yU* are 

the direction cosines (with respect to the vertical) of the photon 
path before and after encounter with the dust particle* Since 
the dust may absorb the radiation as well as redirect it, the true 
scattering probability is 

A — 

P (cos e) = Ci>g 2 (,cos &) 

where the single scattering albedo cOg describes the probability 
that a photon is scattered rather than absorbed in a single 
dust particle encounter ( Wo = 1 for pare — thus conservative — 
scattering) . 

If the scattering particles are spherical and the incident 

A 

light is unpolarized, the phase function P(cos 9j can be computed 

from Hie theory if the complex refractive index m = m^ ~ i mj^ 

and the ratio of the particle radius to the wavelength of the 

incident radiation are known for the particles. When the dust 

layer contains particles of several different sizes, the phase 

function must be weighted by the number of particles having that 

size and then averaged over all possible sizes. Thus, fhe -phase 

function for a homogeneous dust layer where n(r) is the number of 

particles per unit volume whose radii lie between r and r + Sr is 

ORIGINAL PAGE IS 
OP POOR QUALITY 
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given by 


I * 

? (cos s') = T 4 ntr’) cl^, 


CD 


where r, and are the limits of the size distribution and N is 
the total number of particles per unit volume. The dust layer is 
homogeneous if the refractive index and n(r) do not change 
horizontally or vertically throughout the dust layer. Although 
only single scattering has been mentioned above, the optical 
parameters required by the multiple scattering approximations to 
the radiative transfer equations discussed in section III can be 
•generated from the phase function given by Eq. 1. For instance, 
the single scattering albedo for- the size distribution n(r) is 
given by 


So = i / ?((ios &) cUiose) . (2) 

For the calculations presented here we have assumed that the 
particles in the dusty Martian atmosphere are spherical and that 
they can be described by the standard gamma (or two-parameter) 
distribution used by Pang et al . (197 6) ; 


a.b-1 itlb _ X- 


n(r) = N • (ab) *> ' r ** e' / fL '"V*"! 
where I' is the gamma function. Since 

a= gj nfr ncr) dt* 

where G is the geometric cross-sectional area of particles per 
unit volume and is given by 


(3) 


C3a) 

C3b) 


f TT JU-r) dr 
’'f. 


C3c) 
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the parameters a,b can be interpreted as the effective radius and 
the effective variance of the size distribution about a, respectively. 

By varying a,b and the components rij. , of the refractive index. 
Pang ^_t a^, (1976) estimated— on the basis of a least-error curve 

fitting of the 'Mie theory using Eqs. 1 and 3 to the phase functions 
constructed from the 1971 Mariner 9 dust storm data for two 10 nm 
bands centered at 268 and 305 na — that 

0,8 4a 4 1,8 
0.2 < b , . 

By studying these and other bands in the ultraviolet (UV) spectrum 
in the same manner. Pang and Ajello (1976) have determined ’’best® 
values for the refractive indices in the UV. Except for three 
uncertain points with wavelengths less than 0.2 / jm ., their data 
Kavfr been reproduced by the crosses (+) plotted in Pig, 1 at the 
appropriate band centers. 

To extend these results into the visible region, we have used 
the information provided by the isotropic single scattering albedos 
(3^) estimated from the Mariner 9 TV observations of limb haze 
during the Martian global dust storm (E, Anderson, communicated 
by C, B, Leovy) . The values of for the Martian tropics are 

shown in Table I. For optical depths characteristic of the active 
or quasi-steady Martian 1971 global dust storm C 1-2) , these 

£5^ depend only slightly on the total optical depth T‘« . 

A 1*^ ^ 

For isotropic scattering P(cos fe) - (a constant). Thus, 

the phase function given by Eq. 1 and its corresponding are 
independent of particle size. If the actual scattering by the dust 
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is not isotropic, however, the estimated from the data by 
^ssn^i^^g isotropy are size dependent. Assuming that the van de Hulst 
(1974) or two-stream similarity relations are valid, this apparent 
value of the Isotiropic is related to the actual (non-iso tropic) 
single scattering albedo by 



/ -r. 


(4) 


where g is the phase function asymmetry factor defined by 

I = i ?Ccas e),Cos 6 elCcos o) . (5) 

This dependence of the apparent value on parameters generated 
from the actual phase function introduces the size dependence of 



Thus, we proceed as follows: Utilizing the size distributions 

established by Pang et al . (1976), we use Mie theory with Eqs 1 and 
2 to compute Sj at a specific wavelength but for a range of re- 
fractive indices. Then is calculated from Eq. 4. Figure 2 

shows various curves obtained for \ = 0.586yu?n.. Taking a value 

for 3^ from Table I, we find 2^ on a curve generated for the 
desired size distribution and real part of the refractive index. 
From the abscissa of the graph we read the value for m^ and by 
moving vertically to the corresponding curve for , we find the 
actual value for the single scattering albedo. 

In this way we have established the six points indicated by (o) 
in Fig. 1. Because the curves in Fig. 2 have fairly shallow slopes 
with respect to mj_ , the uncertainties tabulated in Table I-yield 
greater uncertainties in m^ (shown by the er'ror bars in Fig. 1), 



7 


Unfortunately, the values listed in Table I depend on an 

absolute interpretation using vidicon data; as such, they are 
not very reliable, (The violet filter data are particularly 
dubious.) In spite of this, these are consistent with values 
estimated by independently analyzing the Mariner 9 B-camera data 
(Leovy 6^ aJ:.* > 1972). Furthermore, they are consistent with 
values obtained from the geometric albedo of the Martian surface 
(primarily bright areas) from which the dust ‘presumably came. 

To compute these latter , we relate to the geometric 


albedo Ag for a homogeneous, isotropic semi-infinite medium: 


I 




I ji.Aji.jXi AjLc 


( 6 ) 


where 


with 

©0 = COs’^ . 


^ " XX (7) 

being the solar intensity arriving at zenith angle 
The diffuse intensity 



£a 

TT 






< 0 * 


( 8 ) 


emerges at an angle 6 = cos’^ytt . The H functions are tabulated 
in Chandrasekhar (I960, p. 139). The resulting curve for Ag is 
shown in Fig. 3. Using de Vaucouleur’s (1964) data for the 
geometric albedo, we find from Fig. 3 and then use graphs 
similar to Fig. 2 to obtain the m^ values marked (a) in Fig. 1. 

We have used the geometric albedo here instead of the Bond albedo 
because Ag is the observed (from the earth) quantity. 

Most of the calculations in this paper were done for one or 
both of the two size distributions listed in Table II, From Fig. 2 
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we see that these two size distributions more or less envelope 
the estimates at visible and near-infrared wavelengths constructed 
from size distributions consistent with' the work of Pang et al . (1976), 

'Similar m^ values have been used for both size distributions. In 
the UV the real part of the refractive index is taken from Pang and 
Ajello’s C1976) work; elsewhere mr is taken to be consistent with 
M,ead ’ s C1970) estimates. The values used are shown in Fig., 1, The 
computed or m^ values at visible and near-infrared wavelengths 
are insensitive to variations of mt^ vithin the range 1,55 im^ 4 1.80. 

To demonstrate the consistency of the various data sources and 
methods, we have plotted cd^ , which is one of the most critical 
optical parameters required to describe the radiative transfer in 
a dusty atmosphere, in Figs. 4 and 5 as constructed for size 
distributions S-II and S-I, respectively. Also plotted are the 3o 
computed from the refractive indices measured for two basalt 
samples, one studied by Pollack al , (1973) and the other by 

Egan et al. (1975). Again, the S-II and S-I size distributions 
were used. Basalt is interest because Adams and McCord (1969) 
have matched the geometric albedo of Mars in the visible and near- 
infrared spectrum with features of the reflectance curves of 
limonite-stained basalt. Differences between the optical properties 
of the two samples shown here are attributed to dlf'ferences in 
chemical composition. 

Figure 1 shows our best estimate for the complex refractive 
index at solar wavelengths. The UV data are taken as they 'stand 
except for the point at X = 0.36^Yit whose low m-^ value was required 
by Pang and Ajello (1976) to correctly model the Martian opposition 

effect. They refer to Mead (1970), but Mead’s work does not by 
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itself impose such a stringently low value. The Mariner 9 TV 
data points are taken to be correct despite their possible bias. 

In the region 0.65 ~ — 1 . , we have used the refractive index 

estimated from our interpretation of the planetary geometric 
albedos. Por I.O^tto, we have required m'^ = 0.003, which is 
suggested by the basalt data. Csee Appendix B.) 

The two different curves in Fig. 1 for Wj, at the visible and 
near-infrared wavelengths are. a direct result of the size dependence 
of mj^ when computed using values derived from spacecraft or 
earth-based observations by assuming that scattering is isotropic 
in the dusty Martian atmosphere. In the following sections 
reference to the size distributions S— I and S— II will generally 
mean not only the distribution’ given by Eq. 3 with the size 
parameters listed in Table II, but will also imply the associated 
ffij. and m^ values from Fig. 1. 
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III. Radiative Transfer 

To compute the vertical distribution of solar heating', we 
shall use simplified but reasonably accurate -models of Che multiple 
scattering radiative transfer equations. Because these models are 
computationally fast, they are exceptionally useful for parametric 
studies and can be easily incorporated into -more general radiative 
equilibrium studies. Furthermore, our lack of deta-iled knowledge 
about the optical properties of the Martian dust and its spatial 
extent necessitates models which require only a few basic 
parameters . 

In Appendix A we consider four simplified multiple scattering 
approximations: a f our-*s tream discrete ordinate model (Liou, 1974), 

a four-stream multiple average-intensity model, and two versions 
of a delta-Eddington model. The accuracy of these are compared 
against an essentially exact doubling method fOr conditions 
appropriate to dusty Martian atmospheres. 

The delta-Eddington -method (Joseph et al . , 1976) is a special 
combination of the usual Eddington and forward peak truncation 
approximations. More explicitly, the delta-Eddington azimuthally 
independent phase function is given by 

'p^gCtose) = + (/-f) ( / + 3 W-S &) (9) 

where' Q is the scattering angle between the incident and emergent 
directions identified by the direction cosines and t respectively, 
while f is the fractional scattering partitioned to the forward 
peak described by the delta function The factor is 
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related to the phase function asymiaetry factor by 



t/: 


Cos e 


?jj(cos e) dCcosft*) - f + (/.;) 


5 ^' 


( 10 ) 


If the optical depth and single scattering albedo are scaled as 


Tj « (/-3 b Or (11) 

« (/-f) S. /{ /-3.I) J (12) 


the delta-Eddington radiative transfer solution can be obtained 
from the standard Eddington approximation simply by replacing 
T* (Ofl » g by , gg , all of ’which can be computed once f is 

known. Since 

f = J e ) P^^Ccos©) d (cos 0) (13) 

where Pjj denotes the ith Legendre polynomial, we can set 
f = if is the coefficient of Pj^ in a Legendre polynomial 

expansion of the normalized phase function. If the actual phase 
function has the same second moment as the Heny ey-Gr eenstein 
function, then f = g*" , 

The results shown in this and subsequent sections have been 
computed using the delta— Eddington method with f “ where g and 
Wa are determined from a Mie calculation using the refractive 
indices of Pig. 1 with the appropriate sise distribution listed 
in Table II. Of the four methods examined in Appendix A» the 
f *= version of the delta-Eddington had the most couaiscontly 
small errors when computing radiative fluxes and flux divergences. 
Purthermore, its overall error was snallest for the 
small solar zenith angles. Each of the other «ethode bad region, 
in Yfji. space where they were as good or better then the 
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Eddingtoa approxination , but the discrete and multiple intensity 
methods in particular tended to have uncomfortably large errors 
for nearly vertically incident insolation. 

In the calculations that follow, the Martian atmosphere was 
assumed to be a plane~parallel atmosphere containing a vertically 
and horizontally homogeneous dust layer extending to the ground. 

The .boundary conditions specified were that the diffuse downward 
radiation vanish at the top of the dust layer ('T = 0) , assuming 
that absorption by in the solar spectral region is negligible, 

and that 

(14) 

where T,; is the total vertical extinction optical depth of the 
dust cloud, Aj is the surface albedo, and is the cosine of 
the solar zenith ang’le, while F4 and FI- are the upward and 
downward radiative fluxes, respectively. 

To determine the surface albedo, we compute the Bond 
albedo for a homogeneous, isotropic, semi-infinite atmosphere: 

f -C' X' 

where S and are given by Eq . 7 and Eq. 8, respectively. 

The resulting curves for and the ratio Ag/A^, known as the 

phase integral, are shown in Fig. 3. The Bond albedo is the 
appropriate surface albedo because it represents the r'atio of the 
flux reflected in all directions from the planar surface to the 
flux intercepted by that surface. Using Fig. 3 and the computed 
for the S-II size distribution (including the associated refractive 
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indices from Fig. 1), we obtain the surface albedo shown in Fig, 6, 
except that we have held the ground albedo constant at Aj = 0.480 
for wavelengths beyond 1.4^irt.. 

Obviously, there are, large uncertainties inherent in this 
derivation for A^. The consistency of the values computed from 
the geometric albedos and of those derived from the Mariner 9 T7 
data suggest that the dust in the Martian atmosphere is the same 
as the dust on the surface. If they do have the same size 
distribution and the same refractive indices, there is no need 
for a surface albedo at all. The low’ere boundary condition given 
by Eq. 14 could be replaced by 

Ft(T,i,^o) = 0* at ca . (16) 

As we shall show in the next section, inserting the Bond albedo 
constructed using Eq, 15 into Eq, 14 with a finite does indeed 
yield practically the same radiative fluxes and heating rates for 
‘T’ti as does utilizing Eq. 16 directly for a much larger , 
However, the dust on the surface may not be the same as the 
dust particles that become airborne. They may have different size 
distributions or the dust on the surface may be mixed with ice. 
Furthermore, the airborne dust in a global dust storm need not 
have originated at the surface immediately beneath it. 

If the two dust regimes are optically different, than Eq. 14 
gives the appropriate lower boundary condition, but the ground 
albedo in such a case may well be different from that shown in 
Fig. 6. 


ORIGINAL 
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IV. Solar Heating 

To examine the vertical distribution of solar heating due 
to dust in the Martian atisosphere, we first define the following 
quantities: 


( Aj ) 

= / If'" dh 

(17) 




<^2 

= •- * Ot 

(18) 


- j a. 

(19) 


where is the vertical extinction optical depth for dust 
evaluated at X = while As represents the interval 

0,2 4 X •£ 5.Q^« and are the heating rates per unit 

volume and per unit mass, respectively. mr radiative 

flux and includes both the direct and diffuse solar contributions, 
while j is the mass density of air (CO^ in this case). Q ,j. is a 
useful measure o-f the total solar heating rate because it is 
independent of any vertical redistribution of the dust as long as 
the total optical depth remains the same. Local heating rates 
(Q 2 ,Qjjj), of course, will vary according to the local dust 
concentration. 

The integral in Eq, 17 was evaluated by a 32*-point Gauss- 
Legendre quadrature for .tn(.2) 4 (| ==>tn.X) A .^.tl (5 , ) . Since 
the global dust storms originate during the Martian southern 
.hemispheric summer, the solar constant was evaluated for a 
planetary distance of 1.45 A.TJ., where it has the value 
64,6 mw cm~^. (1 mw = 10^ ergs sec~^.) The wavelength distri- 

bution of the solar Intensity was taken from Allen(1973) . The 
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optical parameters ( Sa , g, ) were computed for the quadrature 
wavelengths using the refractive indices in’Fig. 1; these parameters 
are listed in Appendix B. Unless otherwise noted, the surface 
albedo AS-11, shown in Fig. 6, was used for all calculations. 

Qc^(Aj) is plotted in Fig, 7 as a function of optical depth 
T and solar zenith angle cos’"^y{e • To insure that the solar 
heating rate being averaged over wavelength refers to a single 
height level, the optical depth for a single wavelength has been 
sca-led with respect to T , the optical depth at X = 0.586yiw. 

f 

Since 


dr 

dt 


<r; 


( 20 ) 


where njj is the number density of dust particles and is the 

extinction cross-section for X= 0,58 ^m. (assumed constant), this 
required that the integral in Eq. 17 be evaluated by summing over 




the quadrature points using values of ^ ^ ^ which 


n - 




(21) 


Thus , 


dr 


^3* _ ^«T d Fritx 

" ~ 


( 22 ) 


The solid curves in Fig, 7 assume that the dust on the 
surface is optically the same as the dust in the atmosphere so 
that the atmosphere-ground system can be treated as a homogeneous 
semi-infinite medium. The other curves show that placing a 
physical surface with the ground albedo shown in Fig. 6 at finite 
optical depths typical of the Martian global dust storm* 

C 1t' 2, Leovy et al , , 1972) changes the solar heating per 
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unit optical depth hy less than 5% for jx^ = 1.0 and by even 

smaller margins for smaller * In fact, for tO.6, the 

curves for = 2,0, 1.5, and 100. agree within 1%. The computed 

net fluxes also change by only 1-5%, and the largest differences 

again occur when yKo — 1.0. This good agreement is not unfounded 

since we have used the same basic assumption— that the surface 

and airborne dust particles are optically the same— -to construct 

the ground albedo curve in Fig. 6. As long as this assumption 

is true, the values on the solid curves will differ from the 

< 

curves computed for a given finite by less than 6% for T i , 

As shown by Fig. 8, tbe relative errors of the infinite 
versus finite T(j curves for the solar heating contributed by 
different segments of the solar spectrum are comparable to those 
for the heating integrated over all of Thus, the agreement 

between the curves shown in Fig. 7 is not a result of the 
cancellation of much larger changes contributed by different 
wavelength intervals. Figure 8 also shows that the largest 
contribution by far to the total solar heating rate comes from 
the visible region C Ayts * 0 *^ 4X4 followed by the 

near-infrared 1.0 A \ t 2.5 ^ir). The region beyond 

contributes very little to the integrated solar heating. 

Comparing curves in Fig. 9 reveals that there is little 
difference between Q,^(Aj) computed from the S-I rather than the 
S-II size distribution even though the contributions from the 
various spectral regions differ significantly. The increased 
visible region contribution for S-I--due to the larger 
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indices associated with, the S-I size distribution at visible 
wavelengths — is almost exactly compensated by the decreased 
near-infrared contribution — due to the generally smaller ratio 
of particle radius to wavelength for the S-I distribution-’-as 
compared to the S-II results. Relative differences in Q.j-(dj) 
computed for these two distributions are less than 3% while the 
corresponding fluxes differ by a similar amount. Usually, 
differences resulting from acceptable changes in the size 
distributions are much smaller than differences arising from 
possible surface albedo changes. As shown in Fig. 10, for instance, 
reducing the ground albedo by 50% dcreases the values for 

/lo = 1.0 by 10% at the top of the dust cloud and by almost 20% 
a t *r “ 1.5. 

because the visible region SO dominates the integrated 
solar heating, replacing the surface albedo AS-II by the constant 
value A„ = 0,30 does not significantly alter the results (within 
4%) for Q^(Ag). The corresponding net fluxes, however, change 
systematically by 5-12%, As one might expect, most of this 
difference is associated with the upward flux . Indeed, 

F'KTf^) decreases by less than 4% when A = .30 replaces AS-II 

3 

% 

while F^(0) decreasesby as much as 15%. The constant surface 
albedo increases the UV contribution to Q^(Aj), but the decrease 
of the near-IR contribution compensates this increase almost 
exactly. 
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We can define a solar heating rate for the entire dust layer 

by ^ 

A. r " 

Q^Cd.'^ cir = • (23) 

The variation of with is shown in Fig. 11. Again, we see 

A 

that the results are quite similar 'when. = 1.5 for values 
computed from either the S-I or S-II size distributions, but that 
is substantially reduced for large over a less -reflecting 
surface. Also shown in Fig, 11 is the layer heating rate when 
is reduced by a factor of 10, 

As shown in "Fig. 12, Q^(dj) varies linearly with r for 

when = 0.15. Although the Q^(As> rates are comparable 
with the rates obtained for Tn - 1.5, the volume and mass heating 
rates will decrease by an order of magnitude because there is less 
dust to do the heating when TJ; = 0.15. 

Figure 12 shows the solar heating rates as a function of 
T and when f a 0.15. Again, the two different boundary conditions 
• given by Eqs. 14 and 16 give similar results. The largest error 
between the two is Figure 13 shows the effect of replacing 

the AS-II ground albedo by the constant A^ = ,30, The relative 
errors have the same sign as for the comparison when = 1»5, 
but are smaller in magnitude by a„ factor of two or so. Reducing 
the ground albedo by half, however, changes the solar heating 
rates by <^20% throughout the dust column. 

To translate the solar heating rates per unit optical depth 
into heating rates per unit volume (Qg) or unit mass (Qjjj) » 
need to know the vertical distribution of the airborne dust. 



19 


Dust was detected as high as 70 km during the 1971 Martian global 
dust storm (AjelLo et al . , 1976). Assuming that the dust was 
uniformly mixed by volume throughout the atmosphere, Eq. 20 can 
be integrated to yield 

'T(z) = c • p 

where p denotes pressure and the constant c is given by 


C * 




?s 


(24) 


where q = njj/n^ is the 'volume mixing ratio, is Martian gravity, 
™C 02 ~ ^ 10“^^ gm per molecule (CO 2 ) , and p^ is the Martian 

surface pressure (taken to be 5 mb). Since the total optical 
depth during the dust storm is 1-2 (Leovy e t al . , 1972), 

c^i>0.3 mb^^. This implies that q 1.3 x 10”"^^ dust particles 
per air molecule for the S-II size distribution and q ~ .8 x 10“'^^ 
for S— I, The number' of dust particles in an atmospheric column 
is then / 9 ^ , which is 1,5 x 10^ cm”^ for the S-I size 
distribution and 2,4 x 10^ for S-II, 

To compute the solar heating per unit mass shown in Eig. 14, 
we multiply the Ty = 1.5 curves in Fig. 7 by ■oaPT' IS 

OP POOR 


• 1 tir 

*'7 31 = » 


which is constant for the the uniformly mixed dust. Thus, the 
solar heating per unit mass is fairly constant with height for 
ju,et 1.0, The decrease near the surface is more pronounced at 
all .as increases beyond til, 5 (see Fig. 7). If we 
denote averages over (global averages) by an overbar (Qjjj, for 
example), we see that the averaged solar heating rate decreases 
significantly in the lower part of the dust layer. Furthermore, 
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Qjjj(^b= 0,4) is a good approximation to the averaged heating rate 
above z = 1.0 (r~0,6) whenT^^ 1.5. 

Specific values of the solar heating per unit mass, Qjji(/ia)» 
and the averaged value are tabulated in Table III for 
representative heights and for = 1*5, 0.15. To convert the 
units given in the table (mw gm”^) to K day”^, the ‘numbers should 
be multiplied by ‘'-100 . Thus, may exceed 50K. day"^ while 
= 1.0) may reach 90K day~^ when = 1.5. Even in the 
optically thin case when T;< = 0.15, heating rates of nearly 
10 k . day”^ may occur when the sun is overhead. 

In Table IV we have shown the partitioning of the incident 
solar radiation into the layer heating (per unit mass), the 
upward flux out of the layer, Ff(0), and the net flux into the 
surface . During the Martian global dust storm some 20% of 

the incoming solar radiation was absorbed by the dusty atmosphere, 
and the flux to the surface was decreased by one-third as compared 
to the optically thin case, (ihis latter result assumes that the 
ground albedo has not been changed by the dynamically active dust 
storm^ The upward flux out of the dusty atmosphere remains the 
same. Thus, we would expect the atmosphere to become warmer 
during a Martian global dust storm, as the Mariner 9 observations 
have shown that it did. We cannot theoretically determine how 
much warmer the dusty atmosphere should become, however, until 
the dust contribution to the infrared cooling is known. 
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V. S^mlmary and Discussion 

Due to radiative scattering, the solar heating per unit 
mass o£ a dusty Martian atmosphere in which the dust is uni- 
formly mixed is surprisingly constant with height if the 
total extinction optical depth ^ 1.5. This result was 
suggested long ago by the Mariner 9 IRIS temperature obser- 
vations of the 1971 global dust storm. The- heating rates 
per unit mass derived in the previous section are as large 
as 90 k/ day when = 1.5, 

The solar heating derived by Moriyama (1975) for the 
Martian global dust storm is considerably larger than that 
computed here and the heating attains a prominent maximum in 
the middle layers of the dust cloud. Moriyama chose m^ = 

0,033 in the visible region since this value produced a single 
scattering albedo of 0.7 for \ 0.6pm. The value 0,7, how- 

ever, assumes isotropic scattering (Leovy et , 1972) and 
the corresponding non-isotropic is close to 0.9 (see Fig. 

2), The increased absorption resulting from this unrealis- 
tically large m. value accounts for his enhanced heating 
maximum in the middle layers of the dust cloud and the de- 
creased solar heating near the ground. Our results do support 
his conclusion, however, that even optically thin dust layers 
can produce heating rates of a few degrees per day. 
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If the van de Hulst similarity relations are valid, the 
largest uncertainty affecting the results given here is 
probably the ground albedo in the visible spectrum, although 
assuming vertical homogeneity may turn out to be equally 
damaging. 

To complete the description of the radiation field in a 
dusty Martian atmosphere, the vertical distribution of the 
infrared cooling by the combined dust-C 02 system needs to be 
calculated. As pointed out by Moriyaraa (1974), most of the 
large solar heating of the dust-laden atmosphere is probably 
balanced by enhanced infrared cooling due to the airborne 
dust. 
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APPENDIX A 

Comparing approximations to the radiative, transfer .equation 

for conditions appropriate to Mars 


We consider here four approximations to the radiative transfer 
equation; a four-stream discrete ordinate method, a four component 
multiple average- intensity method, and two versions of the delta- 
Eddington method described in section III. 

Both delta-Eddington methods, DE(g^) and DEC^Vr)* assume that 
the fraction f of the total scattering partitioned to the forward 
peak of the delta-Eddington phase function (Eq. 7) is related to 
the second moment of the normalized phase function: 





(A() 


Por the DE( ) version, f is computed from the Legendre polynomial 
expansion of the normalized phase function generated by a Mie 
program for the appropriate size distribution and refractive index 
at the given wavelength. The DE(g^) version, however, assumes that 
the actual phase function P(cos 9) has the same second moment as 
the Henyey-Gr eenstein phase distribution: 

f = ^ (H-G) = 

O 


where g is the asymmetry factor defined by Eq, 5 * For both models 
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the asymmetry factor is computed from Mie theory, again using a 
given size distribution and the appropriate refractive index for 
the specified wavelengths. 

In the four-stream discrete ordinate method (denoted FS) the 
angular intensity distribution for a given wavelength, (r,^) , 
is replaced by 2J=4 intensity streams, which are the 1^ evaluated 
at the 2J Gauss-Leg endre quadrature points. Thus, the radiative 
equation of transfer for solar wavelengths becomes discretized 
with respect to the scattered and incident direction cosines 



where is the normalized phase function, W* is the single 

scattering albedo, is the cosine of the solar zenith angle, and 
Ij is the solar intensity. The source function integral has been 
replaced by a 2J Gaussian quadrature with quadrature weights aj at 
-the quadrature points jx.^ . The above equation holds for a specific 
wavelength although the X subscript has been suppressed. When 
J = 1 or 2, algebraic formulas can be obtained for the different 
intensity streams i< = l,2J.(Liou, 1974; note, however, 

that his formula assumes a different lower boundary condition than 
Eq. 1-H in section III.) 

Essentially, the four-stream method assumes that the phase 
function is given by the first four terms of its Legendre polynomial 
expansion : 

a3-fa3 j 

?(cos e) 5‘ X 'PU.u.'') 5. %. (as) 

^ ^ 1U* 
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This severe truncation of the phase function expansion leads to 
negative values of P(cos G) for some scattering angles. This 
unrealistic result may account for the suprisingly large errors 
of the FS model for nearly vertical insolation (yzj"-!.). 

This severe truncation of the phase function expansion for 
a small number of intensity streams can be circumvented if we 
assume that the angular intensity distribution is constant over 
the 2j intervals covering the range -1, £* y«. i +1, We 

can still use the discrete ordinate formulas if we make the 
following substitutions; 









where the quadrature intensity streams are replaced by intensity 
components (*r) , assumed to be the constant value of I(T,yi) on 

the interval ± jj, t. . Since energy conservation requires that 

aJ 

i't 

the multiple intensity method used here assumed that ai , the 

Gaussian quadrature weights* Since oij « -1. and = Py “ 'bbe 

remaining are specified once the are known; this also 

predetermines the median points . 

As shown by Fig, A— 1 through A~4, this multiple average- 
intensity method (denoted FMS) models the shape of the net flux 
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divergence as It varies with optical depth quite well, but the 

nagnitude is too large when compared with the Doubling method 

results. Generally, the FMS error is comparable to that produced 

by the discrete ordinate method. 

Figures A-1 through A-10 compare either fluxes or net flux 

divergences (with respect to optical depth) as computed by the 

radiative approximations with the same quantities computed by the 

Doubling method. The Doubling method used h'ere utilizes 22 terms 

of the Legendre polynomial expansion of the phase function 

generated by a Mie calculation using the S-II size distribution 

for Tjj = 1.5 and the S-I distribution with its slightly smaller 

particles for = 0,15. (In this section the vertical optical 

depth T applies to the given wavelength.) Interchanging the size 

distributions for these two cases does not significantly change the 

results given here as long as the corresponding refractive indices 

are taken from Fig. 1. All calculations shown in the figures, 

except for Fig. A-8, are computed for the wavelength. X= 0.586^m- 

and a surface albedo A = 0.2788, which are representative of the 

s * 

important visible wavelengths. In all cases the incident solar 
flux has arbitrarily been set equal to SF = 1.0, ' 

The two delta-Eddington methods geneially agree quite well 
with each other; in Fig. A-9j for ’instance, they are indistinguish- 
able for J-La - 0.6, Both approximations .are usually significantly 
better than either of the four-stream m'odels — especially for the 
larger values. For ‘T^'^l.O, the intensity stream models are 
distinctly better when computing flux divergences only for 
^ 4 , 0.6, As noted above, both are considerably worse than 
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DE(g^) when ~ 1.0, Figures A-6 and A-7 show relative errors (as 

compared to the Doubling method) of the various approximations when 

computing the downward flux to the surface or when calculating the 

upward flux at the top of the dusty atmosphere. Again, the delta- 

\ 

Eddington methods are generally better for large JS.q , and DECg*^) 

Is significantly closer to the Doubling method result than DE(^^/s) 
for ^a'^1.0. The FS model, however, is quite good for all when 
computing the downward flux at the surface. 

The relative errors when computing the layer heating rate, 
defined as shown in Fig. A-5 for % = 1.5 

and in Fig. A-10 for = 0,15. Both Figs. A-5 and A-7 show the 
nature of the four— stream error, which can become uncomfortably 
large for 

Figure A-8 shows’ the flux divergence computed at )v = O.lOl^w. 

for %i = 1.5 and A = 0.03. Again the delta-Eddingt on approximations 

s 

are best. However, the multiple average-intensity method is 

noticeably better than the discrete ordinate method for the 

stronger absorption (d5o'^0.6) at the ultrayiolet wavelengths. 

The FMS method is far worse than the others for the almost 

optically thin case where T’n = 0.15. When yio “ 1.0, relative 

errors in the flux divergence computed from the four-stream multiple 

. ' 

average-intensity method can reach 35%. The corresponding error 
for the layer heating rate can exceed 30%, as shown In Fig. A-10. 

The four— stream discrete ordinate method, however, is competitive 
with the delta-Eddington approximation when Tjj “ 0,15. 
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Thus, we conclude that for conditions appropriate to the 
dusty Martian atmosphere the delta-Eddington DE(g^) method is 
adequate and is the best of the four methods examined above. This 
is clearly true for the fluxes and flux divergences (with respect 
to optical depth) computed for > 0.6. Furthermore, the DE(g^) 
method requires the fewest number 'of optical parameters (I3e , g, 

, Tjj , and 0^^,^ , the extinction cross-section) of the four 
methods, and it is almost an order of magnitude faster computa- 
tionally than the four-stream methods. 

The four-stream discrete ordinate method is noticeably 
better for while the four-stream multiple average- 

intensity model may be useful for certain long slant optical 
paths or at strongly -absorbing ( tSa i. ,5) wavelengths. 

We again emphasize that the comparisons are for the size 
distributions and refractive indices (at solar wavelengths) 
thought to be appropriate for Martian dusty atmospheres. 

However, the general trends discussed above indicate that the 
delta-Eddington method should be considered whenever one needs 
a computationally fast, simple yet realistic approximation to 
the full radiative transfer equations applied to scattering 

\ 

atmospheres . 
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Appendix B 

Table B-I lists the optical parameters used by the delta- 
Eddington approximation to the radiative transfer equations. 
These parameters were generated using the S~I and S-II size 
distributions from Table II with the tabulated refractive 
indices at the listed quadrature wavelengths. ' The refractive 
Indices listed below (and shown in Fig, 1) were obtained by 
linear Interpolation over wavelength on the values given in 
Table B-II. 
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TABLE B-I 


Optical parameters and constants at the quadrature wavelengths. 

is the product of the quadrature weight and the solar intensity 
at the quadrature wavelength. The solar insolation as a function 
of wavelength was taken from Allen (1973) and adjusted to a 
planetary distance (for Mars southern summer) of 1.45 A.U. The 
32 Gauss-Legendre quadrature points were taken for the interval 
2ytm) A ^ 5 .Oyttw-) . The scaling factor r = 

also listed. All parameters are associated with or generated 
from either the S~I or S-II size distributions. 


S - I 


X(^m) 

fV' 

L3q 

§ 


r 


,r 



.201 

.618 

.845 

.895 

.901 

4.95 

1. 

80 

.02200 

.205 

. 618 

.8 4,4 

.894 

.900 

16.64 

1. 

80 

.02200 

.212 

.620 

,838 

.894 

.900 

41.51 

1. 

81 

.02192 

.222 

.631 

.824 

.898 

.904 

86.42 

1. 

85 

.02126 

.236 

. 644 

.815 

.901 

.907 

141.07 

1. 

85 

.01967 

.255 

.656 

.805 

,908 

.914 

354.73 

1. 

85 

,01877 

.278 

. 690 

.777 

.916 

.922 

920.57 

1. 

85 

.01529 

■ .308 

.741 

.752 

.917 

.923 

3007.70 

1. 

84 

.01134 

.344 

.858 

.715 

.921 

,927 

5614.71 

1. 

75 

.00507 

.388 

.856 

.701 

.947 

.954 

7633.26 

1. 

75 

.00585 

.442 

.855 

.685 

.967 

.974 

16527.19 

1. 

75 

.00630 

.508 

.872 

.675 

,992 

.999 

20789.93 

1. 

75 

.00630 

.586 

.886 

.658 

.993 

1.000 

23149.82 

1. 

75 

.00619 

.680 

.948 

.621 

1.031 

1.038 

23135.47 

1. 

75 

.00273 

.793 

.950 

.580 

1.043 

1.050 

21569.70 

1. 

75 

.00259 

.925 

.959 

.626 

1.056 

1.063 

18816.43 

1, 

65 

.00300 

1.081 

.965 

.620 

1.082 

1.090 

15888.35 

1. 

65 

.00300 

1.262 

.969 

,611 

1.101 

i:i09 

13137.76 

1, 

.65 

.00300 

1.470 

.973 

,610 

1.089 

1-097 

10650.49 





1.706 

.977 

.605 

1.089 

1.097 

7936.04 





1.970 

.979 

. 607 

1.055 

1.062 

5213.59 





2.261 

.982 

.607 

1.007 

1.014 

' 3791,23 





2.575 

.983 

. 607 

.944 

.'951 

2359.92 





2.907 

.984 

,605 

.871 

.877 

1578.33 





3.249 

.985 

.602 

.792 

.798 

1136,04 





3.592 

.985 

.597 

.715 

.720 

816,64 




. 

3.924 

.986 

.591 

.643 

. 648 

502,18 





4.232 

,986 

,585 

,581 

.585 

344.50 

> 

r 


r 
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TABLE B-I (cont.) 


4.503 

.986 

.578 

.530 

.534 

244.42 

1.65 

.00300 

k.llh 

.986 

.573 

.492 

.495 

157.84 

1.65 

.00300 

4.886 

.986 

.569 

.465 

.468 

89.12 

1.65 

.00300 

4.978 

.986 

,566 

.451 

.454 

34.90 

1.65 

.00300 


S - II 


.201 

.579 

.892 

5.76 

.919 

4.95 

1 . 

80 

.02200 

.205 

.580 

. 890 * 

5.76 

.919 

16.64 

1 . 

80 

.02200 

.212 

.582 

.888 

5.77 

.920 

41,51 

1 . 

81 

.02192 

.222 

.589 

.879 

5.79 

.923 

86.42 

1 . 

85 

.02126 

.236 

.597 

.873 

5.81 

.926 

141.07 

1 . 

85 

.01967 

.255 

.606 

.866 

5.83 

.930 

354.73 

1 , 

85 

.01877 

.278 

.634 

.845 

5.86 

• .935 

920.57 

1 . 

85 

,01529 

,308 

.681 

.816 

5.88 

.938 

3007.70 

1 . 

84 

.01134 

.344 

.812 

.770 

5.94 

.948 

5614.71 

1 . 

75 

.00493 

.388 

,847 

.749 

5.98 

.953 

7633.26 

1 . 

75 

.00415 

.442 

.871 

.730 

6.05 

.966 

16527,19 

1 . 

75 

.00367 

.508 

.885 

.715 

6.18 

.986 

20789.93 

1 . 

75 

.00359 

.586 

.899 

.701 

6.27 

1.000 

23149.82 

1 . 

75 

.00344 

.680 

.955 

.660 

6.39 

1.019 

23135.47 

1 . 

75 

.00157 

.793 

.958 

.634 

6.61 

1.054 

21569,70 

1 . 

75 

.00149 

.925 

.960 

.628 

6.62 

1.056 

18816.43 

1 , 

75 

.00188 

1.081 

.955 

.618 

6.88 

1.097 

15888.35 

1 , 

73 

.00225 

1.262 

.956 

, 611 

7.12 

1.135 

13137.76 

1 . 

70 

.00257 

1.470 

.962 

. 614 

7.29 

‘ X ,163 

10650.49 

1 . 

66 

.00295 

1.706 

.966 

. 606 

7.66 

1:222 

7936.04 

1 . 

65 

.00300 

1.970 

.971 

.604 

7.85 

1 . 25’3 

5213.59 

1 . 

,65 

.00300 

2.261 

.976 

.605 

8.00 

1.277 

3791.23 

1 . 

,65 

.003 

iOO 

2.575 

.979 

.609 

8.05 

1.285 

2359.92 





2.907 

.981 

.612 

7.98 

1.273 

1578.33 





3.249 

,983 

.615 

7.76 

1.238 

1136.04 





3.592 

.984 

. 616 

7.44 

1.187 

816,64 





3.924 

.985 

.616 

7.05 

1.125 

502.18 





4.232 

.986 

.613 

6.66 

1.062 

344.50 





4.503 

.986 

.611 

6.30 

1.004 

244.42 





4.724 

.986 

. 608 

' 5.99 

.956 

157.84 


/ 



4.886 

.986 

.606 

5.78 

.921 

89.12 


f 



4.978 

.986 

.604 

5.65 

.901 

34.90 

1 , 

.65 

.00300 
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TABLE B-II 


Refractive indices for the solar wavelengths 0.2^w. -X - 5,0y«w< 
were generated from the following values using linear interpolation 


over wavelength X 

. The 

TJV data (first section) are 

from Pang 


and Ajello 

(1976) 

. All 

values 

listed here are shown 

in Fig 

. 1. 


X C^m. ) ! 


i.210 

.220 

,233 .243 .253 .269 

.280 

.305 

.326 

Hj. (S— I and 

S-II) ; 

1.80 

1,85 

1.85 1.85 1.85 1.85 

1.85 

1.85 

1.80 

m^CS-I and 

S-II) : 

.022 

.0215 

,020 .019 ,019 .017 

.015 

.012 

.007 

X (ynw)i 

.340 

.414 

.585 

,610 .700 .800 

.900 

1.00 

il.50 


1.75 

1.75 

1,75 

1.75 1.75 1.75 

1.65 

1.65 

1.65 


1.75 

1.75 

1.75 

1,75 1.75 1,75 

1.7'5 

1.75 

1.65 

: 

.0050 

.0063 

.0063 

.0039 .0024 .0026 

.0030 

.0030 

.0030 


.0050 

.0037 

.0035 

.0022 .0014 .0015 

- 

.0021 

.0030 
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Table I 

l.sotropic single scattering albedo u3^ estimated from the 
Mariner 9 TV observations of limb haze over the Martian 
f-Topics. It has been assimied that 1. (E. Anderson, 

fcommunicated by C, Bo Leovy.) 


B-Camera Filter 

Effective Wavelength 

• 


Violet 

. 414pm 

0o62 

± .05 

Orange (1) 

,585pm 

0.72 

o 

• 

-H 

Orange • (2) 

. 610pm 

0.81 

± .05 
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TABLE II 

Parameters used in Eq. 3 to generate size distributions for the 
solar heating calculations. These parameters are within ranges 
established by Pang e_t al. (197 6), The geometrical cross-sectional 
area per particle, calculated from Eq. 3c, is also listed. 


S-I : 

a = 1*0 

b = 0.4 

6 = 0.376 


S-II: 

a X • 3 

b = 0.25 

G = 2.649 
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TABLE III 


Solar heating rates per unit mass (mw gm"^) for a uniformly 
mixed dusty atmosphere. 


S-II, 

AS-II, i; = 1.5 


Cos 

(solar 

zenith 

angle) 

- 


Average 

p(mh) 

z(scale ht.) *t(X 

= . 58 6^m) 

0.2 

0.4 

0.6 

0.8 

0.9 

1.0 

over 

.333 

2.7 

0.1 

.52 

. .66 

.74 

.80 

.82 

.85 

.58 

2.0 

0.92 

0.6 

.22 

-47 

.64 

-77 

.83 

.88 

: .45 

5.0 

0. 

1.5 

.10 

.29 

.49 

.68 

.77 

,86 

,36 

S-I, 

AS-II, r^= 0.15 




.333 

2.7 

0.01 

.061 

.068 

.074 

.080 

.083 

.086 

.063 

2.0 

0.92 

0.06 

.054 

.065 

.073 

.080 

.084 

.087 

.059 

5.0 

0. 

0.15 

.044 

.060 

.071 

.080 

.085 

.089 

.054 
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TABLE IV 


Radiative fluxes, shown as per cent of the incident solar radiation, 
for a uniformly mixed dusty atmosphere. 


S-II,AS-II, 1.5 

Cos 

0.2 

(solar zenith angle) = 
0.4 0.6 0.8 0.9 

1.0 

Average 
over /it 


25% 

24% 

21% 

19% 

19% 

18% 

21% 

Ft(0) 

52% 

45% 

40% 

36% 

34% 

33% 

39% 


23% 

31% 

39%- 

45% 

47% 

49% 

40% 

5^- Ig (mw cm~2) 

12.9 

25.8 

38.8 

51.7 

58.1 

64.6 

32.3 

S-I, AS-II, -Tjj =0.15 

5b = 

6% 

3% 

2% 

2% 

2% 

2% 

3% 

Ft (0) 

43% 

40% 

38% 

37% 

37% 

37% 

38% 


51% 

57% 

60% 

61% 

61% 

61% 

59% 
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Figure 1 


Figure 2 


List of Figures 

Real and imaginary parts of the complex re- 
fractive index m = - im^. The points (+) 

were taken from Pang and Ajello (1976), The 
points (o) are for a basalt sample studied by 
Egan ^ al, (1975); the points (x) are for a 
similar sample studied by Pollack £l, (1973), 
The remaining m^ points were evaluated using the 
two different size distributions listed in Table 
IL. The points (^) are based on the Mariner 9 
TV data while the points (□) are based on the 
planetary geometric albedo. The values are 
consistent with Mead’s (1970) work. 

Generating the isotropic and nqn- isotropic 
single scattering albedos for X = .586 as a 
fianction of the imaginary component of the re- 
fractive index. Four size distributions were 
used: solid lines are for an effective radius 

a = l.Opm, dashed lines for a = 1,5pm, while 
(o) indicates an effective variance b = 0.4, 
Otherwise, b = 0.25. (The parameters a, b were 
used in the formula given by Eq. (3) with N = 1 
particle per unit vol\ima.) = A65", 
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Figure 3 


Figure 4 


Figure 5 


Figure 6 


Geometric and Bond albedos, and Aj, respectively, 
as a functio'n of the isotropic single scattering 
albedo in a homogeneous, isotropic, semi-infinite 
medium. The phase integral IP = is also 

shown. 

Single scattering albedo 6i)o generated from different 
sources by various methods using the S-II size 
distribution. 

Single scattering albedo generated from different 
sources by various methods using the S-I size 
distribution. 

The ground albedo constructed assuming that the 
airborne dust is optically the same as the dust on 
the surface. The circles (o) are for the ground 
albedo estimated (1*7 linear interpolation) at the 
quadrature points used for the integration over 
the solar spectrum. Collectively these are designated 
AS-II. The upper curves show 6Ja at the quadrature 
points as calculated for the S-1 and S-II size 
•distributions using the refractive indices shown in 
Fig. 1, These points have been joined by straight 
line segments to show the general trend. 
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Figure 7 


Figure 8 


Solar heating per unit optical depth as a function 
of T, the optical depth at \ = .586^m, and 
the cosine of the solar zenith angle. Q (a ) 

T S 

has been computed for S-II optical parameters and 

various combinations of surface albedo A and 

s 

total optical depth The units are 1 mw = 

io4 -1 

10 ergs sec . Integration over wavelength 
was accomplished by using a 32 point Gauss -Legendre 
quadrature over In \ . 

Contributions to the total solar heating (integrated 
over A ) from the following regions: 

Auv : 0-2 1 X ^ O.SVni 

' ^NUV' ^ ^ 0.4^m 

Ay^g* ^ X ^ 1,0pm 

Ajjis= 1-0 1 X 1 2.5pm 
Ag^ 0.2 ^ X ^ 5.0pm 

The effects of the lower boundary conditions given 
by Eqs . 14 and 16 are also shown. The heating rates 
were computed for p^ = 1.0 and S-II optical parameters. 
The contributions fpr the various regions were 
obtained by summing the values for those quadrature 
points falling in the specified region a. . 

A 
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Figure 9 


Figure 10 


Figure 11 


Figure 12 


Figure 13 


Contributions to the total solar heating rates 
from the same spectral regions as defined in 
Fig, 8, Calculations were done for = 1.0, 

Tj^ = 1.5, and for both the S-I and S-II optical 
parameter sets . The surface albedo used was 
taken from Fig . 6 . 

The effect of surface albedo on the total solar 
heating rate for various values as a function 
of optical depth. Calculations were done for 
Tjjj = 1.5 and for S-II optical parameters - 
The layer heating rate as a function of 
zenith angle for both t>, = 1.5 and t„ = 0.15. 

The effects of various optical parameter sets 
and surface albedos are shown. 

Same as Fig. 7, except that ~ 0.15 is the . 
’’finite" optical depth and the S-I parameter 
set was used. 

Comparing the contributions to the total heating 

from various spectral regions when = 1.0 and 

= 0.15. The effect of replacing the albedo 

from Fig. 6 by a constant = 0.3 value is shown. 

Note that when = 0.3, the contribution 

is indistinguishable from the A contribution. 

NIR 
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Figure 14 


Figure A-1 


Solar heating rates per unit mass as a function 
of height for an atmosphere in which dust is 
uniformly mixed and 1.5 (at = 5 mb). 

The calculation used the S~II optical parameters 
and AS-II ground albedo. The broken line shows 
the value of the solar heating rate (per unit 
mass) averaged over 

Comparing different approximations to the radiative 
transfer equation with the Doubling method for 
X = .586|jm, A^ = .2788, and = 1.0. The plotted 
quantity is the divergence of the net flux mth 
respect to the optical depth at \ = .586pna. The 
solar flux (SF) has been arbitrarily set equal to 
one. The total optical depth is = ^*5, and 
the optical parameters were computed from the S-II 
size distribution. 

• Doubling method 

DE (f = g^) 

DE (f = 6 ) 2 / 5 ) 

_ FS (4-stream discrete ordinate method) 

FMS (4-stream multiple average -intensity 

method) 
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Figure A-2 
Figure A-S 
Figure A-4 
Figure A-5 


Figure A- 6 


Figure A- 7 


Figure A- 8 


Figure A-9 


Same as Fig, A-l except = 0.9. 

Same as Fig. A-1 except = 0.6. 

Same as Fig, A-1 except = 0.2. 

Relative differences (in °L) of the radiative transfer 
approximations -when compared against the Doubling 
method for computing the heating rate for the entire 
layer, defined as 

% = I = ^NET<^N> - 

A ® 

Q is shown as a function of u » the cosine of 

T ro 

the solar zenith angle . The calculations were 

done for - 1.5, X ~ .586)jm, A^ - ,2788, and 

the S-II set of optical constants. 

Same as Fig. A-5, except the quantity being plotted 
is the net flux at the surface (Tj^ = 1.5). 

Same as Fig. A-5, except the quantity being plotted 
is the upward flux Ft (0) from the dust layer. 

Same as Fig. A-1 except plotted for A^ = 0,3 and 
X = ,201^m. 

Comparing different radiative transfer approximations 

with the Doubling method for dFj^g.j/dx generated 

using the S-I optical constants, = 0.15 and 

A = 0.2788 for \ = . 586pm. 
s 

Same as Fig. A-5, except computed using the 
S-I optical constants and for* =* 0.15. 


Figure A-10 
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